Entanglement dynamics in chaotic systems 
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We study quantum chaos for systems with more than one degree of freedom, for which we present 
an analysis of the dynamics of entanglement. Our analysis explains the main features of entan- 
glement dynamics and identifies entanglement-based signatures of quantum chaos. We discuss en- 
tanglement dynamics for a feasible experiment involving an atom in a magneto-optical trap and 
compare the results with entanglement dynamics for the well-studied quantum kicked top. 
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Theories of chaos and of quantum mechanics juxtapose 
in the discipline of 'quantum chaos' (QC), which has at- 
tracted significant theoretical and experimental 0, 0- 
efforts. Recently the advent of quantum information |4| 
has highlighted the role of entanglement as a resource, 
and stimulated theoretical studies of entanglement in QC 
systems [!, IH IE El E3 • Here we present a general anal- 
ysis of entanglement dynamics for unitarily evolving QC 
systems, and apply this analysis to a system of significant 
experimental interest: an atom in a magneto-optical lat- 
tice (AMOL), which we show is a feasible experimental 
system to observe and test entanglement dynamics and 
to rapidly enhance entanglement for certain initial states. 
The theoretical methods and results are general, which 
we demonstrate by application to the well-studied quan- 
tum kicked top (QKT) @, |H H El ■ 

We also identify 

the initial entanglement rate as a signature of quantum 
chaos. 

Entanglement features of QC systems are rather sub- 
tle: in fact recent theoretical studies reveal a dichotomy: 
entanglement may be enhanced as the chaoticity parame- 
ter is increased |fj,Q or the opposite can occur y|. In fact, 
this dichotomy can be understood by observing that the 
former result applies for a global chaotic quantum system 
for which entanglement applies between its constituents, 
whereas in the latter case one considers weak coupling 
between a quantum chaotic system and the environment 
which could be another quantum chaotic system. Our fo- 
cus is on the former case which can yield a rapid increase 
of entanglement. 

In particular, we analyze the evolution of an 
AMOL [J^l under realistic experimental conditions and 
show that entanglement arises between the atomic spin 
and motional degrees of freedom. The ability to tomo- 
graphically reconstruct the reduced density matrix of the 
atomic spin state [l5j makes it possible to study the evo- 
lution of this entanglement experimentally. We analyze 
the dynamics of entanglement (via the linear entropy of 
the spin density matrix) and characterize the global dy- 



namics by the size of the chaoticity parameter and local 
dynamics by whether the initial state is supported pri- 
marily by regular or chaotic eigenstates of the Hamilto- 
nian. AMOL experiments would allow the first empirical 
studies of entanglement evolution in a QC system. 

An AMOL provides an attractive framework for study- 
ing entanglement evolution in a QC system with and 
without coupling to an environment. This is because, 
in addition to having more than one degree of freedom 
and the ability to tomographically reconstruct states, de- 
coherence can be controlled by detuning the laser relative 
to the atomic resonance frequency. Decoherence is negli- 
gible for the far off-resonance AMOL so coupling to the 
environment can be ignored (unitary dynamics), and en- 
tanglement between the spin and motional subsystems 
of the overall chaotic system can be explored. By tuning 
the laser frequency close to atomic resonance, coupling 
to the environment is increased and the resulting effects 
of entanglement can be observed. 

Here we are concerned with the former case of the far 
off-resonance magneto-optical lattice for which coupling 
to the environment can be neglected and entanglement 
is enhanced. We show that under feasible conditions, 
the AMOL can exhibit generic features of entanglement 
dynamics, for example quasiperiodicity for a state ini- 
tially localized in a regular regime and a rapid increase 
of entanglement with no subsequent quasiperiodicity in 
a chaotic regime. 

Generic features of entanglement in QC systems can 
be understood by examining the spectral properties 
of the evolution operator U(t) on the system Hilbert 
space H = ®f =1 Ti {i) with H {i) the Hilbert space of 
dimension di for the i th subsystem. We consider two 
common categories of unitary evolution: (i) U(t) = 
exp(—iHt/h) for a time- independent Hamiltonian H and 
(ii) U[t = tit) = F n with F being a Floquet oper- 
ator F = Texp[— i/h J Q H(t)dt] and t a discrete time 
variable. The unitary evolution operator can be ex- 
pressed as U(t) = Y Ji exp(-iu)it)\^i)((/)i\ for {\4> t )} a 
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time-independent orthonormal basis of Ji of dimension 
d and {exp(— the corresponding eigenvalues. The 
evolution of an arbitrary initial density operator p(0) 
over time t is U(t)p(0)W{t) = ^ e" 1 ^ l^><<^l for 
Wij- = w< - and ee (0^(0)1^-). 

Entanglement for pure states (p = p 2 ) with two sub- 
systems is given by the entropy of the reduced density 
operator p(t) of either subsystem. The linear entropy 
5=1— Tr(p 2 ) is a convenient measure of entanglement, 
with S = for no entanglement and 5 = 1 — 1/di for 
maximum entanglement. The time-dependent entropy is 

S(t) = 1 - ]T C iiW e- i ^« +w «) t (1) 

i,j,k,l 

with 

Qjkl = PijPkl ( U P> v m\(l)i)(4>j\Up,V n ) 

m,n,p,q 

(Ug,V n \(j)k)(<pl\u q ,V m ) (2) 

for {|w,„)} orthonormal bases for and . 

Diagonal elements pa quantify support of p on [/- 
eigenstates . These eigenstates can be associated with 
regular and chaotic regimes 0, 0] ; hence a state p can 
be identified with classical regular or chaotic regimes (or 
a combination) by its support on [/-eigenstates. We em- 
ploy this correspondence between support 0] and reg- 
ular vs chaotic dynamics to characterize entanglement 
dynamics for an AMOL and other QC systems. 

The AMOL is subjected to a far off-resonance AC 
Stark shift as a function of atomic position z given by [lij ] 

4 

V = —V\ cos0l cos2fcz — fi ■ B c s (z) , (3) 
o 

for fi = fvyF = —psF/F, 7 the gyromagnetic ratio and 
F the total angular momentum vector of the hyperfine 
ground state. As an example, we consider 133 Cs with 
F = 4 and psB c s(z) — ~V\ sin Ol sin 2kz e z + B x e x for 
k the laser wave vector, O^ the relative polarization an- 
gle between the counter-propagating laser beams forming 
the 1-D lattice and V\ the single-beam light shift. The 
coupling between the spin precession and center-of-mass 
motion leads to entangled spinor wave packets and B x is 
the tunable chaoticity parameter. (In fact, Eq. (4) is ap- 
plicable to more general systems if the periodic potential 
is replaced by a harmonic potential 0.) In the classical 
limit, Eq. (3) describes a magnetic moment interacting 
with the same B e g(z) 0, with F/F the direction vec- 
tor for the classical angular momentum, and the classical 
four-dimensional phase space parametrized by atomic po- 
sition and momentum (z,p) and direction (9, </>), of F/F. 

We seek to study entanglement dynamics for lattice 
parameters that are accessible in current experiments. 
Therefore we choose V\ = 160-Er, Ol — 80° and [IbB x = 
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FIG. 1: Classical Poincare sections for the AMOL at E — 
~280E R for Vi = 1G0E R , Q l = 80° and fi B B x = 12E R with 
(a) p y = 0, d^iy/dt > and (b) p = 0, dp/di > 0. 
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FIG. 2: Entanglement S versus time r = Ent/h in the AMOL 
for an initial state that is localized on (a) a regular island with 
(z/X,p/hk,9,cj>) = (-0.15,0,1.27,0) and (b) the chaotic sea 
with (z/\,p/hk,9,4>) = (0.06,0,7r/2,0). The inset shows the 
initial increase of entanglement for the regular (solid) and 
chaotic (dashed) initial states. 

12E R , for E R = H 2 k 2 /2M, the recoil energy. Classical 
Poincare sections for these parameters and total energy 
E = p 2 /2M + V = —280E R reveal a mixed phase space 
with islands of regular motion embedded in the chaotic 
sea (Fig. 1). Quantum states are localized to phase space 
coordinates (z,p, 9, <ft) by preparation in a product of the 
motional and spin coherent states \a = z + ip}\9, (j>) |2(lj . 

An AMOL state localized around (z, 0, 9, <fi) can be pre- 
pared by cooling atoms to the ground state of the diabatic 
potentials. The lattice is then shifted until this state is 
centered at (z/X). The spin is rotated until the Bloch 
vector is pointing in the direction (9, (f>). We pick an ini- 
tial state that is centered on an island in the Poincare 
section in Fig. 1. For comparison we also consider an 
initial state in the chaotic sea. 

The evolution of entanglement, quantified by linear en- 
tropy, for states that are initially regular (Fig. 2(a)) or 
chaotic (Fig. 2(b)), exhibit two main signatures of chaos. 
At short times, entanglement in the chaotic regime in- 
creases at a faster rate than for the regular regime, 
thereby supporting the concept that chaos can cause 
rapid generation of entanglement as observed in other 
systems 0, 0|. Also oscillations are prevalent for ini- 
tially regular states but not for chaotic states (which has 
been also observed for the quantum kicked top 8]). We 
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FIG. 3: Population pa vs corresponding eigenenergy Ei = fkOi 
for the AMOL with an initial state localized in (a) a regular 
island with (z/\,p/hk,6,<p) = (-0.15,0,1.27,0) and (b) the 
chaotic sea with (z/\,p/hk, 6, <j>) = (0.06, 0, tt/2, 0). 
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FIG. 4: (a) For the initial state localized on the regular island 
the entanglement dynamics (solid curve) can be reproduced 
by only considering the evolution of the four main pairs of 
eigenstates (dashed curve). The long term behavior (b) shows 
quasiperiodic motion with multiple frequencies. 



explain how the power spectrum of S can provide a sig- 
nature of chaos, as Lahiri suggested |9j, by exploiting the 
(7-eigenbasis. 

U -eigenstate support for the states initially in the reg- 
ular and chaotic regimes are depicted in Fig. 3. The ini- 
tial state on the regular island has support dominated by 
four pairs of regular eigenstates with each pair nearly de- 
generate. This support over few eigenstates is responsible 
for the quasiperiodic evolution of linear entropy, and Fig. 
4(a) shows that excellent replication of entanglement dy- 
namics is possible by only including these four pairs of 
eigenstates. In contrast, the initial state in the chaotic 
sea has support over a larger number of the chaotic set of 
eigenstates extending over a broader frequency spectrum 
(Fig. 3(b)), due to a breakdown of semiclassical theory 
in the chaotic regime 0, 0] . 

A rigorous understanding of entanglement evolution 
emerges by noting that S(t) in Eq. (JJJ depends on 
eigenfrequency-difference sums u>ij + LUki , which can be 
identified in the power spectrum of S(t). The fast os- 
cillations in Fig. 2(a) are due to the large differences 8u 
between the four main peaks in Fig. 3(a), and small fre- 
quency differences between almost degenerate eigenstate 
pairs at each peak in Fig. 3(a) result in slow oscillations 
with long periods as seen in the long-term behavior of 
the entanglement in Fig. 4(b). The terms + lj^i that 
appear in the evolution correspond not just to differences 
in the eigenfrequencies, but also can be a sum of w,j +lum- 
For example, the main oscillation in Fig. 2(a) is due to 
the sum of two difference frequencies. 

The key feature of dynamical entanglement for our pur- 
pose is the initial increase of entanglement for a chaotic 
state, which is more rapid than for the regular state (inset 
of Fig. 2). This rapid rise in entanglement is not unique 
to the AMOL system 0, Q . The rate of increase of en- 
tanglement is obtained by expanding Eq. (1) at t = 0, 
which reveals a quadratic increase as a function of time, 
S = (t/t ) 2 with t — 0.01; this behavior is surprising 
at first because an exponential increase is expected for 
states in the chaotic regime and a quadratic increase for 
states in the regular regime |l0|. Of course this expec- 



tation applies for the asymptotic semiclassical regime, 
but our system lies in a deeply quantum regime and 
hence the chaotic state is not well localized, with non- 
negligible support over the regular regime which results 
in the quadratic increase in entanglement. For the initial 
state in the chaotic regime, Ax/X « 0.07, Ap/hk w 2.7, 
Api x = and Afj, z — Afi y = l/\/2. 

Our methods to analyze the AMOL are generic and 
generalize to other unitarily evolving QC systems as we 
now show for the quantum kicked top with Hamilto- 
nian [HEDS 

OC 

H=^-Jl+pJ v ty-nr), (4) 

J n— — oo 

for J x ,J y ,J z su(2) operators and k the chaoticity pa- 
rameter. The QKT can be constructed from a collec- 
tion of N — 2j qubits in the symmetric representa- 
tion with collective spin operators J a = 53i=i ^t 21 ; ano - 
{(Jia) the Pauli operators for the i qubit 8]. For 
k = 3, t = l,p = tt/2 entanglement behavior is similar to 
that of the trapped atoms described here || . Bipartite 
entanglement between a pair of qubits and the remaining 
qubits reveals quasiperiodic evolution for an initial state 
centered on an elliptic fixed point. For a state centered in 
the chaotic sea, no quasiperiodic motion is present, and 
just as in our AMOL, an initial rapid increase of entan- 
glement is observed confirming this generic behavior. 

Support of an initial state over (U — F")-eigenstates, 
of the Floquet operator F = exp(— iAtJ^/2jr)exp(— ipJ y ) 
is shown for an initial state centered on the elliptic fixed 
point and one in the chaotic region of Ref . ||| . The initial 
state centered in the regular region can be mainly decom- 
posed into a few (size y/N) 'regular' eigenstates 0. In 
this case the state which we have localized at a fixed point 
has most of its support on three eigenstates, of which two 
are degenerate. The difference between the correspond- 
ing eigenphases <f> m , 5<ft = 0.003 determines the frequen- 
cies of oscillation in the evolution of S (Fig. 6(a)). 

The flat entanglement power spectrum for the chaotic 
state is due to the broad support (size TV) of the 
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FIG. 5: Population p mm vs corresponding eigenphases <f> m for 
the QKT with an initial state localized in (a) a regular island 
and (b) the chaotic sea in Fig. 1 of 0. 



crease of entanglement for initial states in the chaotic 
regime. Our analysis relies on studying support over the 
unitary evolution eigenbasis, which applies generically to 
other unitarily evolving QC systems, as we demonstrate 
with the quantum kicked top. Thus we have introduced 
a means for extending QC experiments to more than one 
degree of freedom, observing and understanding entan- 
glement dynamics in such systems, exploiting quantum 
chaos for a rapid increase in entanglement, and resolving 
the dichotomy between cases where chaos enhances vs 
diminishes entanglement generation for initially chaotic 
states. 
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FIG. 6: (a) For the regular initial state of the QKT the entan- 
glement dynamics (solid curve) can be reproduced by consid- 
ering the evolution of the three highest weighted eigenstates 
(dashed curve) in Fig. 5(a). (b) The rise time for the state in 
the chaotic sea is roughly exponential. 



initial state on the 'chaotic' [/-eigenstates (Fig. 5(b)). 
Because N — 50 qubits is in the semiclassical regime, 
the distinction between regular and chaotic entanglement 
dynamics is more pronounced than what we observed for 
the AMOL, which was not as semiclassical. (As we con- 
sider 133 Cs for the AMOL, F = 4 is fixed and thus not 
in the semiclassical regime). Moreover the rise time for 
the initial chaotic state is exponential \l0\ which can be 
regarded as a signature of quantum chaos. 

The QKT system of N qubits behaves collectively like 
a system with one degree of freedom unlike the AMOL 
which has two degrees of freedom: spin and motion. Cou- 
pling between two QKTs, however, allows for entangle- 
ment dynamics between two coupled degrees of freedom 
to be observed, but entanglement is suppressed in the 
strong chaotic regime Q ■ This reduction of entanglement 
can be explained by noting that our results imply that 
the entanglement between qubits in one QKT is chaos- 
enhanced, which causes the qubits to be collectively less 
entangled with the qubits of the other QKT. 

In conclusion, we have presented entanglement dynam- 
ics for an experimentally feasible QC system of atoms 
trapped in a magneto-optical lattice. For realistic exper- 
imental parameters, quantum signatures of chaos exist in 
the dynamics of entanglement, specifically in the initial 
rise and the power spectrum, even when the system is 
not in a semi-classical regime, and there is a rapid in- 
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